Quantum transport thermometry for electrons in graphene 
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We propose a method of measuring the electron temperature T e in mesoscopic conductors and 
demonstrate experimentally its applicability to micron-size graphene devices in the linear-response 
regime (T e w T, the bath temperature). The method can be especially useful in case of overheating, 
T e > T. It is based on analysis of the correlation function of mesoscopic conductance fluctua- 
tions. Although the fluctuation amplitude strongly depends on the details of electron scattering in 
graphene, we show that T e extracted from the correlation function is insensitive to these details. 

PACS numbers: 73.23.-b, 72.15.Rn, 73.43. Qt, 81.05.Uw 



Graphene is an atomically thin graphite layer [l|, Q 
recently used in field-effect transistors Q- In graphene- 
based semiconductor devices phonons are poorly coupled 
to the environment since the mass of carbon atoms is 
typically smaller than that of atoms in the underlying 
substrate, making the overheating of graphene structures 
a likely event at high currents. This raises a question of 
how to measure the temperature of electrons in graphene. 
Since classical conductivity in graphene has a very weak 
temperature dependence at low and intermediate tem- 
peratures [H, extracting the electron temperature from 
transport measurements requires analyzing more subtle 
quantum effects. One possibility would be to analyze the 
decoherence rate t^ 1 using the weak-localization (WL) 
effects in magneto-resistance @, However, it was 
shown theoretically 0| and confirmed experimentally Q 
that the WL in graphene reveals itself in a rather compli- 
cated way due to the influence of inter-valley scattering 
and the disorder which breaks the sublattice symmetry. 
Thus WL does not offer an easy way of measuring the 
electron temperature T e . Another possibility would be to 
exploit the temperature dependence of the amplitude of 
universal conductance fluctuations (UCF) [H, 9, 101. Un- 
fortunately, a quantitative implementation of such anal- 
ysis is hindered by the necessity to both account for 
the temperature dependence of t^ 1 and attribute a def- 
inite symmetry class to a particular graphene-based de 
vice 



11, 12, [3]. However, it has been noticed that 
the correlation functions of random UCF dependence on 
magnetic field B and the Fermi energy ep provide use- 
ful information about subtle spectral characteristics of a 
disordered conductor 
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FIG. 1: The electronic temperature T e , determined via 
Eq. (J), as a function of the bath temperature, T, for four 
graphene flakes (Fl: open circles, F2: open squares, Bl: 
filled circles, B2: filled squares, see Table [IJ. The inset shows 
A c /2.7 for sample D with Lt ~ L x , i.e. not in the regime ©; 
the dashed horizontal line shows the Thouless energy, /i/td. 



fluctuations, 5G = G — (G), as a function of the Fermi 
energy (using a wide range of the UCF magnetofinger- 
prints for statistically representative averaging [lj|): 
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F(A) = ((SG(e F )6G(e F + A))) . (1) 



In this Letter we propose a method of correlation func- 
tion thermometry of mesoscopic conductors and demon- 
strate experimentally its applicability to graphene-based 
micron-size devices. Although the UCF in graphene dif- 
fer in detail [111, EH, EH from those in other mesoscopic 
conductors, the proposed method is robust and indepen- 
dent of such details. It is based on the analysis of the 
normalized correlation function F n (A) of conductance 



This function depends on the electron temperature T e 
which governs the thermal broadening of the Fermi dis- 
tribution. We show that for a quasi-lD wire the width 
A c of F n (taken at the half-maximum) is defined by T e : 



2.7k B T e 



(2) 



This expression allows one to determine T e by extract- 
ing A c from measuring the correlation function ([1]). The 



result is valid with accuracy of about 10% provided that 
the following conditions are fulfilled: 



L T < min(L y , L x ) , 
L y < min(L y , L x ) , 



(3a) 
(3b) 



where = y/ nJD j 'k B T e is the thermal smearing length 
and L v = ^Dt v is the dephasing length, D is the diffu- 
sion coefficient. Inequality pap defines an experimentally 
relevant "high-T" regime (which may extend well below 
IK) and is paramount for the method to work. Inequality 
(|3bj) is less demanding: one can use Eq. ([2|) to determine 
T e also for a 2D sample with aspect ratio a = EJ L y ~ 1 
up to L v /L y ~ 1, albeit with lesser accuracy |16j |. 

We have experimentally tested the method in four 
graphene flakes in the regime ([3]), using low-current mea- 
surements to avoid electron overheating. Then T e should 
coincide with the bath temperature T. We show the re- 
sults in Fig. [H where T e xtracted from Eq. (J5J) and the 
directly measured T are, indeed, in good agreement. 

Below we first derive our main result, Eq. [[2]). then 
proceed with its numerical testing and finally discuss ex- 
perimental results in more detail. 

The brackets ((•••)) in Eq. ((T|) stand for both the en- 
semble and thermal averaging. This equation can be ex- 
plicitly represented [3, llOj as the following convolution: 



F(A) = J deK(e,A)F(e) 7 
K(e, A) = (jp) J dEf(E, e F )f'(E + e,e F + A) 



(4) 



Here F{e) = (SG(E) 8G(E + e)) is the ensemble-aver- 
aged correlator of conductance fluctuations at different 
energies and K is the thermal broadening factor, where 
the energy derivative of the Fermi distribution function 
is f'(E,e F ) = -I /(Ak B T e ) cosh- 2 [{E - e F )/(2fc B T e )]. 
The standard diagrams for T(e) in the lowest order 
in hj {p F t) <C 1 {I is the electron mean-free path) are 
shown in Fig. [2 Structurally, they coincide with the di- 
agrams describing mesoscopic fluctuations in usual con- 
ductors y, |l?| but the so-called Hikami boxes are differ- 
ent Ill.ll2l.ll3| because of graphene-specific features: the 
linear dispersion law, chirality of the carriers and valley 
degeneracy. These features, being paramount for aquan- 
titative description of the UCF in graphene 11,12, 13| , 
have no impact on calculating F n . We will show this by 
analyzing first a narrow graphene wire, Eq. (|3b[) . with 
a strong inter- valley scattering (induced, e.g., by atom- 
ically sharp disorder). Then the part of the correla- 
tor T which contributes to F n can be written via the 
"valley-singlet" diffusion propagators Vnm , neglecting 
the "valley-triplet" modes (see Refs. @,[ll| for the appro- 
priate classification of the diffusion modes in graphene): 



FIG. 2: The diagrams which contribute to the main order in 
the diagrammatic expansion of the conductivity-conductivity 
correlation function. The wavy lines stand for the propagators 
of the diffusion modes describing motion at length scales 2> i\ 
shaded blocks stand for Hikami boxes, describing motion at 
length scales ~ I. 
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, (5b) 



where a = L x /L y and td = L x /D. Under the condi- 
tions © the sum in Eq. (|S"a)) is dominated by the m = 
term. There we have taken into account only the diffusion 
modes and neglected the Cooperons. This corresponds 
to the regime of suppressed WL, when the magnetic flux 
through the area of order L 2 , is much bigger than the flux 
quantum. In this regime a wider range of statistical data 
is available due to the averaging over magnetic fields [l~5| . 

In order to get an asymptotic analytical expression for 
the correlation function F n (A) we assume that in addi- 
tion to Eq. (|3|) the dephasing is sufficiently strong, i.e. 
L x 3> L v . Keeping only the term with m = and per- 
forming the summation over n, we arrive at 
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, (6) 



where t = yj (er^/h) 2 + 1. The correlator T(e) in Eq. ([6]) 
is a sharply peaked function of e with maximum at e = 
and width h/T v . In contrast, the thermal broadening fac- 
tor K(e, A) in Eq. has a broad peak around e = A 
with the width of the order of k B T e 3> fi/r^, according to 
Eq. (|3a|l . Therefore the integration over e in Eq. ([4]) can 
be performed using the mean value theorem, i.e. taking 
K(0, A) out of the integral. As a result the normalized 
correlation function becomes independent of the micro- 
scopic details contained in T{e): 



K{0 7 A) 3(0coth0-l) 



if (0,0) 



sinh 



2k B T e 



(7) 



The width of this function at the half-maximum is 9 C — 
1.36 which results in Eq. ([2]). We stress again that this 
result is truly universal: a precise form of !F(e) is irrele- 
vant for F n (A) in Eq. ([7|). The only requirement for its 
validity is that the function T(e) in Eq. (U) is sharply 
peaked compared to K(s,A). This remains valid under 
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FIG. 3: The normalized correlation function, F n (A), for a 
wire under conditions J3). The solid line corresponds to the 
asymptotic Eq. (J7J). The dotted, dashed, and dash-dotted 
lines correspond to the numerical integration of Eqs. Q with 
a = 0.01, P = 0.3; a = (3 = 0.1; a = 0.1, /3 = 0.3, respectively 
(a = Lt/L x , (3 = Lt/L v ). The width at half maximum is 
A c w (2.8 ± 0.1) fc B T e . The inset shows F n (A) for a square 
sample where each line has the same values of a and /3. 



the condition (|3a[) for any dephasing, L v < L^, and al- 
lowing for all the diffusion modes in graphene or, indeed, 
in any other mesoscopic disordered conductor. 

We have checked this numerically, calculating F n (A) 
for a wide range of T and t v , 0.01 < Lt/L^ < 0.3 
with L v < L Xl and at various values of the symmetry 
breaking parameters. We have also considered a case of 
smooth disorder, taking into account the valley-triplet 
diffusion channels [neglected in the analytical calcula- 
tions in Eqs. (©-((6j)]. A few representative examples 
are plotted in Fig. [3l The values of A c lie within a nar- 
row interval, 2.7 < A c /fcBT e < 2.9, close to the asymp- 
totic value 2.7 of Eq. This shows that the pro- 
posed method works in the quasi-lD case with accuracy 
of about 10%. We have also performed a similar analysis 
of 2D samples, calculating F n (A) for the same range of 
parameters as for wires. The results plotted in the inset 
to Fig. [3] indicate that A c i=s 3fcBT e , i- e - the method still 
can be applied albeit with accuracy of about 25%. 

We have tested the feasibility of the proposed method 
in graphene-based structures. To this end, we have 
experimentally determined the width of the correla- 
tion function ([T]) and thus the electron temperature T e , 
Eq. ([2]). We have compared T e to the bath temperature 
T in low-current measurements, i.e. in the regime when 
graphene is not overheated, and found them to be in sat- 
isfactory agreement, Fig. [1] The experimental samples 
used for correlation- function thermometry are monolayer 
graphene flakes created by mechanical exfoliation [l8[ on 
a n + Si substrate covered by 300 nm of SiC>2. The flakes 
are connected electrically by two Au/Cr contacts in the 
circuit shown in the inset to Fig. |4] and thermally an- 
chored within the 3 He pot of a sorption-pumped cryostat. 
The ac current driven through the sample was 1 nA (at 





L x (/an) 
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L v (/mi) 
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1.4 


1.4 
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70 
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Fl 


4.1 


1.8 
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70 


1.4 


F2 


3.8 


1.8 


>L X 


120 


0.7 


Bl 


3.7 


0.3 


>L X 


80 


0.9 


B2 


2.0 


0.3 


>L X 


50 


0.9 



TABLE I: Characteristics of the samples used for the correla- 
tion function thermometry, Fig. [1] The values of L v are given 
at T = 0.26K; for the last three samples, where L v > L x , it 
cannot be determined from the WL measurements. 



~ 10 Hz). Heating by the current was detected by mea- 
suring the mesoscopic fluctuations, and this current was 
reduced until it had no effect on the fluctuations. The 
Fermi energy was controlled by a gate voltage V g applied 
between the substrate and the flake [l9j]. Samples char- 
acteristics are described in Table fl] The dephasing rate 
in all samples was determined in a way similar to that 
in 0] from a fit of the magnetic field dependence of the 
sample conductivity to the theory of weak localization in 
graphene 0|. For these samples at high carrier density, 
~ 10 12 cm~ 2 , min(L x , L v ) satisfies the conditions ([3]) for 
the applicability of the method. 

Mesoscopic fluctuations of the conductance G occur 
in all samples as a function of both magnetic field and 
Fermi energy, Fig. Ufa), and are reproducible for the for- 
ward/backward sweeps of £f and B. The amplitude and 
correlation function both depend on the bath tempera- 
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FIG. 4: (a) Typical fingerprint of 5G (normalized by e 2 /h) 
for sample Bl at T = 0.25 K and B = 90 mT as a function 
of Fermi energy (controlled by the gate voltage [19] ■) (b) 
Correlation function for sample Bl at T = 0.25 K. The inset 
sketches the circuit used in the experiment. The resistor has 
resistance much greater than that of the graphene flake to 
maintain a fixed ac current of 1 nA. 
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ture T over the full temperature range of the experiment 
(0.25 — 20 K). A range of Fermi energies is chosen such 
that the average resistance does not change significantly 
over this range and contains sufficient number (> 100) 
of fluctuations for averaging (at T > 10 K such a re- 
quirement cannot be satisfied so we restrict analysis to 
T < 10 K). We determine F n (A), Fig. [1(b), and its width 
at half maximum A c performing measurements over the 
magnetic field range 50 mT < B < 300 mT l 15]. It is 
shown in Fig. [I] as a function of the bath temperature 
and agrees with the theoretical relation, Eq. (J2J), with 
T e w T. For comparison, the inset to Fig. [1] shows that 
A c for sample D (which is not in the "high-temperature" 
regime of Eq. pa])) saturates at H/td- 

In conclusion, we have proposed to use the correlation 
function of mesoscopic fluctuations in disordered samples 
for determining the electron temperature from its width. 
We have shown that the method is universal and inde- 
pendent of microscopic details of a disordered mesoscopic 
sample: Eq. ([2|) holds for any sample in regime ([3a| . de- 
termining T e with good accuracy in quasi- ID samples, 
(|3b[) . We have confirmed the viability of the method per- 
forming measurements on graphene devices in the low- 
current regime, when T e « T. For future applications, 
the method may be especially useful for graphene devices 
at higher currents, where overheating is likely to arise 
from inefficient thermal contact with the environment. 
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EP/D031109, the Lancaster-EPSRC Portfolio Partner- 
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